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ABSTRACT 

We  report  the  results  of  a  globally  coordinated  photometric  campaign  to  search  for  transits  by  the 
P  ~  30  d  and  P  ~  60  d  outer  planets  of  the  3-planet  system  orbiting  the  nearby  M-dwarf  G1  876.  These 
two  planets  experience  strong  mutual  perturbations,  which  necessitate  use  of  a  dynamical  (four-body) 
model  to  compute  transit  ephemerides  for  the  system.  Our  photometric  data  have  been  collected  from 
published  archival  sources,  as  well  as  from  our  photometric  campaigns  that  were  targeted  to  specific 
transit  predictions.  Our  analysis  indicates  that  transits  by  planet  “c”  (P  ~  30  d)  do  not  currently 
occur,  in  concordance  with  the  best-fit  i  =  50°  co-planar  configuration  obtained  by  dynamical  fits 
to  the  most  recent  radial  velocity  data  for  the  system.  Transits  by  planet  “b”  (P  ~  60  d)  are  not 
entirely  ruled  out  by  our  observations,  but  our  data  indicate  that  it  is  very  unlikely  that  they  occur. 

Our  experience  with  the  G1  876  system  suggests  that  a  distributed  ground-based  network  of  small 
telescopes  can  be  used  to  search  for  transits  of  very  low  mass  M-stars  by  terrestrial-sized  planets. 

Subject  headings:  stars:  planetary  systems  -  stars:  individual  (G1  876)  -  planets  and  satellites:  general 


1.  INTRODUCTION 

The  G1  876  planetary  system  ranks  with  the  most 
remarkable  dicoveries  that  have  emerged  from  the  first 
decade  of  extrasolar  planet  detections.  As  described  by 
Rivera  et  al.  (2005),  Laughlin  et  al.  (2005),  Butler  et  al. 
(2001)  Marcy  et  al.  (2001),  Delfosse  et  al.  (1998)  and 
Marcy  et  al.  (1998),  there  are  a  number  of  reasons  why 
the  G1  876  planets  are  extraordinary.  Planet  “d” ,  with  its 
1.974  day  orbit  and  M  sin(i)  =  5.9 is  the  lowest-mass 
exoplanet  known  to  orbit  a  nearby  main-sequence  star 
(Rivera  et  al.  2005).  Planets  “c”  and  “b”  are  the  only 
Jovian-mass  companions  known  to  orbit  an  M-dwarf,  and 
their  clear  participation  in  2:1  mean- motion  resonance 
gives  important  dynamical  clues  to  the  formation  and 
evolution  of  the  system  (e.g.  Lee  and  Peale  2002).  Fur¬ 
thermore,  the  G1  876  planets  are  the  closest  exoplanets 
that  have  been  reliably  characterized.  Their  radial  veloc¬ 
ity  amplitudes  induce  an  overall  signal-to-noise  of  nearly 
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100,  and  they  have  have  been  observed  for  well  over  a 
decade. 

The  detection  of  G1  876  d  by  Rivera  et  al.  (2005) 
underscored  the  rapid  development  of  the  Doppler  tech¬ 
nique.  The  detection  of  planets  having  only  a  few 
Earth  masses  is  the  culmination  of  advances  described 
by  Butler  et  al.  (1996),  Marcy  &  Butler  (1998),  Marcy, 
Cochran,  &  Mayor  (2000),  Marcy  et  al.  (2005),  and  Lo- 
vis  et  al.  (2006).  Nevertheless,  while  radial  velocity  pre¬ 
cision  has  improved  to  better  than  1  m/s,  Keplerian  or¬ 
bital  fits  to  Doppler  velocities  give  estimates  of  Mp\  sin(i) 
rather  than  Mp\.  In  order  to  ascertain  the  true  mass  of 
a  planet,  one  must  obtain  an  independent  measure  of 
the  orbital  inclination,  i,  which  is  most  easily  measured 
when  the  planet  is  observed  to  transit  across  the  face  of 
the  parent  star. 

The  now-celebrated  companion  to  HD  209458  provided 
the  first  example  of  a  transiting  extrasolar  planet,  and 
this  body  has  generated  an  incredible  variety  of  obser¬ 
vational  and  theoretical  investigations.  HD  209458b  was 
initially  detected  with  the  radial  velocity  method,  and 
was  then  discovered  to  transit  on  the  basis  of  follow-up 
photometric  observations  at  the  predicted  transit  times 
(Henry  et  al.  1999,  2000;  Charbonneau  et  al.  2000).  HD 
209458b’s  well-sampled  light  curve,  when  combined  with 
high-precision  radial  velocity  data,  has  permitted  pre¬ 
cise  measurements  of  the  planet’s  basic  attributes  (e.g. 
Mazeh  et  al.  2000,  Brown  et  al.  2001,  Laughlin  et  al. 
2005,  Winn  &  Holman  2005). 

To  date,  nine  additional  planets  have  been  found  to 
transit  their  parent  stars,  and  the  varied  uses  to  which 
these  transits  have  been  put  are  described  in  the  review 
article  by  Charbonneau  et  al.  (2006).  Our  goal  in  this 
paper  is  to  ascertain  whether  the  two  outer  planets  orbit¬ 
ing  G1  876  (IL  Aqr,  GJ  876,  HIP  113020)  can  similarly 
be  observed  in  transit.  A  positive  detection  would  give  a 
further  major  improvement  in  our  understanding  both  of 
this  specific  system  and  of  planetary  properties  in  gen¬ 
eral. 

Rivera  et  al.  (2005)  carried  out  a  photometric  search 
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for  transits  by  planet  d  (P  ~  1.974d),  and  showed  that 
such  transits  do  not  occur.  To  date,  photometric  searches 
for  transits  by  planet  c  (P  ~  30d)  and  b  (P  ~  60d)  have 
not  been  reported.  Although  the  a  priori  geometric  tran¬ 
sit  probabilities  for  G1  876  c  and  b  are  only  Ptr  =  1.4% 
and  Ptr  =  0.9%,  respectively,  the  expected  photometric 
transit  depths  are  greater  than  10%.  Librational  motions 
arising  from  the  resonance  between  the  two  outer  planets 
lead  to  predicted  transit  epochs  for  the  outer  planets  that 
are  not  spaced  evenly  in  time.  (The  resonant  structure  of 
the  system  is  discussed  further  in  Laughlin  &  Chambers 
2001;  Rivera  &  Lissauer  2001;  Kley  2003;  and  Beauge, 
Michtchenko,  &  Ferraz-Mello  2005). 

Benedict  et  al.  (2002)  used  HST  astrometry  to  derive 
an  orbital  inclination  for  G1  876b  of  i  =  84±6°.  Their  re¬ 
sult  thus  suggests  that  planet  b  has  a  mass  equal  or  close 
to  the  RV-determined  value  of  M  sin(i)  =  1.89Mj|jp, 
and  that  the  system  is  within  1-sigma  of  being  viewed 
edge-on.  The  implication  of  the  Benedict  et  al.  (2002) 
measurement  is  that  occurence  of  transits  is  much  more 
likely  than  the  a-priori  geometric  probability  would  sug¬ 
gest.  The  Benedict  et  al.  (2002)  result  is  in  conflict, 
however,  with  the  dynamical  models  of  the  system  pre¬ 
sented  by  Rivera  et  al.  (2005),  who  find  a  best  fit  to 
the  radial  velocity  data  when  the  system  (assumed  co- 
planar)  is  inclined  by  i  ~  50  ±  3°  with  respect  to  the 
plane  of  the  sky.  Therefore,  if  transits  are  detected,  then 
the  Benedict  et  al.  (2002)  results  will  be  confirmed  at 
the  expense  of  the  conclusions  of  Rivera  et  al.  (2005). 

G1  876b  and  G1  876c  are  members  of  an  as- yet  un¬ 
observed  class  of  giant  planets  with  intermediate  sur¬ 
face  temperatures.  The  expected  effective  temperatures 
of  Teffb  =  160k  and  Tesc  =  200K  of  the  planets  are 
much  lower  than  the  value  Teg  ~  1350K  measured  with 
Spitzer  for  HD  209458b  (Charbonneau  et  al.  2005),  but 
still  considerably  higher  than  the  Teff  ~  135K  measured 
for  Jupiter.  Determination  of  their  physical  properties 
would  thus  provide  a  useful  link  between  Jupiter  and  Sat¬ 
urn  on  the  one  hand,  and  objects  such  as  HD  209458b 
on  the  other.  Accurate  measurement  of  their  radii  would 
also  give  clues  to  their  interior  structure  (see  Charbon¬ 
neau  et  al.  2006). 

Finally,  the  G1  876  system  with  its  strongly  gravita¬ 
tionally  interacting  planets  and  large  planet-to-star  ra¬ 
dius  ratios,  has  the  potential  to  display  an  extraordinar¬ 
ily  informative  set  of  transit  light  curves.  Our  model 
light  curves  herein  show  what  one  can  expect  when  in¬ 
teracting  transiting  systems  are  discovered,  for  example, 
with  ongoing  radial  velocity  surveys  or  with  space-based 
missions  such  as  Kepler  (Borucki  et  al.  2003)  or  COROT 
(Baglin  2003),  and  allow  for  a  more  complete  understand¬ 
ing  of  the  multiplanet  dynamics. 

The  plan  of  this  paper  is  as  follows.  In  §2,  we  de¬ 
scribe  our  dynamical  model  of  the  system,  and  illustrate 
some  of  the  photometric  properties  that  the  in-transit 
light-curves  would  be  expected  to  display.  In  §3,  we  de¬ 
scribe  an  evaluation  of  the  Doppler  radial  velocity  mea¬ 
surements  that  have  been  taken  during  predicted  transit 
intervals.  In  §4,  we  evaluate  a  range  of  archival  and  newly 
obtained  photometric  data  sets  for  evidence  of  transits. 
In  §5,  we  discuss  our  results,  as  well  as  the  current  status 
and  outlook  for  the  ongoing  TransitSearch  collaboration, 
and  in  §6,  conclude. 


2.  THE  DYNAMICAL  MODEL 

G1  876  (M4V)  is  the  fortieth-nearest  star  to  Earth, 
with  a  Hipparcos-determined  distance  of  4.69  pc  (Per¬ 
ryman  et  al.  1997);  Tycho-II  and  UCAC  positions  vali¬ 
date  this  determination.  As  discussed  in  Laughlin  et  al. 
(2005),  current  estimates  of  its  mass  and  radius  stand  at 
M*  =  0.32  M0,  and  R*  =  0.3  R0.  We  adopt  these  values 
as  fixed  in  this  paper.  To  date,  155  Doppler  velocities 
measured  at  the  Keck  telescope  and  16  velocities  mea¬ 
sured  at  the  Lick  telescope  have  been  published  for  the 
G1  876  system  (Marcy  et  al.  2001;  Rivera  et  al.  2005). 
The  two  planets  induce  a  total  velocity  half-amplitude 
of  the  star  of  nearly  0.25kms_1.  Furthermore,  as  dis¬ 
cussed  in  Rivera  et  al.  (2005)  the  system  displays  a  very 
low  level  of  stellar  radial  velocity  “jitter”,  which  further 
aids  in  obtaining  a  detailed  characterization  of  the  or¬ 
bits.  The  system  has  now  been  observed  for  more  than 
eighty  orbital  periods  of  the  middle  planet,  c.  This  ex¬ 
tensive  time  baseline  allows  for  a  much  more  detailed 
study  of  planet-planet  interactions  than  can  be  obtained 
with  the  other  exoplanetary  systems  known  to  be  in  2:1 
resonance,  a  list  which  now  includes  HD  82943  (Mayor 
et  al.  2004,  Lee  et  al.  2006),  HD  128311  (Vogt  et  al. 
2005),  and  HD  73526  (Tinney  et  al.  2006).  These  three 
systems  all  have  inner  planet  periods  of  order  ten  times 
longer  than  G1  876  c.  Our  baseline  co-planar  i  =  90° 
orbital  model  was  obtained  by  Rivera  et  al.  (2005)  from 
a  dynamical  fit  to  the  Keck  radial  velocity  data.  For 
reference,  the  parameters  of  this  model  are  reproduced 
in  Table  l.13  Additional  fits  with  i  <  90°  were  also  ob¬ 
tained  to  the  Keck  data  using  the  algorithm  described 
by  Rivera  et  al.  (2005).  As  with  the  i  =  90°  baseline 
model,  the  system  is  assumed  to  contain  three  planets, 
all  orbiting  the  star  in  the  same  plane.  We  varied  the 
inclination  of  the  orbital  plane  to  the  plane  of  the  sky 
from  i  =  90°  to  i  =  89.0°  in  decrements  of  0.1°,  and  fit¬ 
ted  for  each  of  the  13+1  remaining  parameters.  N-body 
integrations  of  the  fit  then  generate  a  sequence  of  transit 
predictions  for  the  outer  two  planets.  Uncertainties  in 
the  predicted  transit  midpoints  were  generated  with  the 
bootstrap  Monte-Carlo  technique  described  by  Rivera  et 
al.  (2005).  For  both  outer  planets,  the  uncertainty  in 
the  transit  midpoint  times  was  found  to  vary  from  ~1 
hr  to  a  few  hours.  The  variation  in  uncertainty  from 
transit  to  transit  is  plotted  in  Figure  1.  The  predicted 
transit  midpoints  are  best  constrained  during  the  epoch 
spanned  by  the  radial  velocity  observations.  Errors  in¬ 
crease  as  the  transit  predictions  are  generated  for  earlier 
or  later  times.  Light  curves  were  produced  with  a  pho¬ 
tometric  model  of  the  transit  which  assumes  that  the 
planets  are  opaque  disks  with  radii  R  =  0.93  Rjup  for 
c  and  R  =  1.04  Rjup  for  b.  These  equatorial  radii  are 
computed  using  the  models  of  Bodenheimer,  Laughlin, 
&  Lin  (2003).  The  specific  intensity  of  the  stellar  disk  is 
modeled  with  a  linear  limb  darkening  law,  given  as 

W/A1)  =  (1) 

where  I  is  the  specific  intensity  of  the  stellar  disk, 
v  =  0.724  is  the  corresponding  limb  darkening  coefficient 

13  Note  that  the  elements  reported  in  Table  1  are  osculating 
orbital  elements  (see,  e.g.,  Murray  &  Dermott  1999)  referenced  to 
the  epoch  JD  2452490. 


3 


appropriate  to  V-band  observations  of  an  M4V  primary 
(per  Claret  2000);  and  /i  =  cos(0). 

Figures  2  and  3  show  a  selection  of  model  light  curves 
for  planets  b  and  c.  For  planet  b,  which  has  an  eccentric 
(e  =  0.22),  rapidly  precessing  orbit  ( dw/dt  ~  44°  yr-1), 
the  shape  of  the  predicted  curves  depends  strongly  on 
the  osculating  periastron  angle.  For  i  <  90°,  the  transits 
are  both  shorter  and  deeper  when  the  transits  occur  near 
periastron.  The  orbital  eccentricity  of  planet  b  is  more 
nearly  circular,  with  e  =  0.025.  In  this  case,  there  is 
little  variation  in  the  predicted  light  curves  as  the  orbits 
precess.  The  known  transiting  planets  have  transit-to- 
transit  intervals  AT)  that  are  evenly  spaced  to  within  the 
current  resolution  of  the  observations.  HD  209458b,  for 
example,  has  had  individual  transit  midpoints  measured 
to  an  accuracy  of  several  seconds  (Brown  et  al.  2001, 
Wittenmyer  et  al.  2005),  and  an  average  period  mea¬ 
sured  to  an  accuracy  of  83  ms  (Wittenmyer  et  al.  2005, 
Winn  &  Holman  2005).  If,  however,  one  is  able  to  mea¬ 
sure  differences  in  the  time  intervals  between  successive 
transits,  one  can  derive  additional  dynamical  informa¬ 
tion.  As  discussed  by  Miralda-Escude  (2002),  Holman 
&  Murray  (2005),  and  Agol  et  al.  (2005),  measurements 
of  small  variations  from  exact  periodicity  can  be  used 
to  infer  the  presence  of  additional  bodies  in  the  system, 
potentially  with  the  masses  in  the  range  appropriate  to 
terrestrial  planets.  In  the  event  that  either  G1  876b  or 
G1  876c  were  to  transit,  one  would  expect  the  departures 
from  strict  periodicity  to  be  large.  In  fact,  successive 
transits  of  planet  c  would  occur  at  intervals  differing  by 
as  much  as  4  hours!  By  fitting  a  dynamical  model  to  such 
transit  intervals,  one  could  obtain  information  regarding 
all  of  the  orbital  elements  of  the  planets  (including  incli¬ 
nations  and  nodes). 

3.  SEEKING  EVIDENCE  OF  TRANSITS  IN  THE  RADIAL 
VELOCITY  DATA  SET 

During  a  planetary  transit,  the  occulting  disk  of  the 
planet  passes  over  portions  of  the  stellar  disk  contain¬ 
ing  a  varying  component  of  line-of-sight  velocity  arising 
from  stellar  rotation;  this  phenomena  produces  an  asym¬ 
metry  in  the  star’s  spectral  lines  known  as  the  Rossiter- 
McLaughlin  effect  (Rossiter  1924,  McLaughlin  1924). 
The  effect  has  been  described  in  detail  in  the  context 
of  planetary  transits  by  Ohta,  Taruya  &  Suto  (2005), 
and  has  been  analyzed  for  HD  209458  by  Queloz  et  al. 
(2000),  Bundy  &  Marcy  (2000),  and  Winn  et  al.  (2005), 
and  for  HD  149026  by  Wolf  et  al.  (2006).  In  fact,  the 
transiting  planet  HD  189733  was  initially  discovered  to 
transit  not  via  photometry,  but  rather  via  such  in-transit 
radial  velocity  variations,  found  by  Bouchy  et  al.  (2005). 
For  our  study  of  G1  876,  we  compared  the  observation 
epochs  of  all  of  the  Lick  (see  Marcy  et  al.  2001)  and 
Keck  (see  Rivera  et  al.  2005)  radial  velocity  observa¬ 
tions  to  see  if  any  of  these  observations  occurred  during 
a  predicted  transit  (assuming  i  =  90°).  For  G1  876,  we 
used  a  rotational  period  that  was  determined  to  be  96.7 
days  using  rotational  modulation  (as  reported  in  Rivera 
et  al.  2005);  we  used  this  value  to  calculate  the  am¬ 
plitude  of  the  Rossitter-McLaughlin  effect.  Four  such 
observations  were  found,  three  taken  during  a  predicted 
transit  by  planet  “c” ,  and  one  taken  during  a  predicted 
transit  by  planet  “b”.  Figures  4  and  5  display  the  rel¬ 
evant  parts  of  the  radial  velocity  model  (both  with  and 


without  the  Rossiter-McLaughlin  effect),  with  the  obser¬ 
vations  shown.  The  radial  velocity  model  curves  in  both 
figures  were  generated  by  adding  template  curves  which 
use  the  stellar  and  planetary  radii  and  limb  darkening 
law  given  above  (and  which  apply  the  analytic  integrals 
given  in  Ohta  et  al.  2005),  to  our  fiducial  i  =  90°  three- 
planet  radial  velocity  model.  For  the  transit  epoch  in 
the  case  of  planet  c  (Figure  4),  the  data  is  inconsistent 
with  a  transit.  In  Figure  5,  for  the  outer  planet  b,  the 
situation  is  unresolved. 

4.  A  SEARCH  FOR  TRANSITS  IN  PHOTOMETRIC  DATA 

SETS 

We  have  also  compared  our  model  light  curves  to  a 
variety  of  photometric  datasets  for  G1  876  taken  since 
1989.  The  results  of  these  comparisons  are  as  follows. 

Hipparcos:  The  Hipparcos  Mission  (see  Perryman  et 
al.  1997)  produced  67  accepted  photometric  measure¬ 
ments  for  G1  876  (HIP  113020).  These  measurements 
have  a  median  magnitude  Hp  =  10.148,  with  a  =  0.0045. 
There  is  evidence  for  a  weak  periodicity  when  the  data 
is  folded  at  the  96.7  d  rotational  period  of  the  star  (also 
noted  in  Rivera  et  al.  2005).  In  Figure  6,  we  plot  our 
model  i  =  89.6°  light  curves  for  G1  876  b  and  c  super¬ 
imposed  on  the  Hipparcos  epoch  photometry.  None  of 
the  Hipparcos  measurements  were  made  near  or  during 
a  predicted  transit  window. 

Fairborn  Observatory: 

A  number  of  long-term  photometric  monitoring  pro¬ 
grams  are  being  carried  out  with  the  automated  tele¬ 
scopes  at  Fairborn  Observatory  (e.g.,  Henry  1999;  Eaton, 
Henry,  &  Fekel  2003)  and,  as  described  in  detail  in  Rivera 
et  al.  (2005),  G1  876  has  been  on  one  of  the  observing 
programs.  Figures  7  and  8  show  the  Fairborn  Observa¬ 
tory  data  and  our  corresponding  model  light  curve  for 
planets  b  and  c.  Data  taken  during  predicted  transit 
epochs  for  planet  c  show  no  indication  of  a  transit.  The 
Fairborn  Observatory  data  indicate  that  long-term  stel¬ 
lar  variability  from  G1  876  is  of  order  0.05  magnitude, 
which  is  somewhat  smaller  than  the  variation  observed 
with  Hipparcos. 

TVansitsearch  and  AAVSO:  Seagroves  et  al.  (2003) 
describe  Transits earch,  a  cooperative  distributed  observ¬ 
ing  project  involving  sub-meter  class  telescopes  world¬ 
wide.  The  TransitSearch  strategy  is  to  observe  known 
planet-bearing  stars  at  the  dates  and  times  when  tran¬ 
sits  are  expected  to  occur.  It  is  therefore  a  targeted 
search,  which  differentiates  it  from  ongoing  wide-held 
surveys.  We  identify  observing  windows  for  candidate 
stars  through  the  use  of  the  bootstrap  Monte  Carlo  tech¬ 
nique  described  by  Laughlin  et  al.  (2005).  The  observa¬ 
tional  campaigns  are  prioritized  by  the  a-priori  likelihood 
that  a  particular  candidate  planet  will  display  transits. 
This  likelihood  is  given  by, 
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where  a  is  the  semi- major  axis  of  the  orbit,  i?*  is  the 
radius  of  the  star,  Rpi  is  the  radius  of  the  planet,  zu  is  the 
longitude  of  periastron,  and  e  is  the  orbital  eccentricity. 

The  first  TransitSearch  dataset  for  G1  876  was  obtained 
by  one  of  us  (Shankland)  from  2003  Oct.  27.9583  to  2003 
Oct.  28.0368,  with  a  V-hltered,  16-bit  speckle  ccd  at  the 
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primary  focus  of  a  40.4  cm  f/5.1  reflector  (at  air  mass 
from  sec(z)  =  2.7  to  1.7).  The  photometric  time  series 
shows  a  monotonic  flux  increase  of  0.16  +/-  0.04  mag 
from  2003  Oct.  27.9583  to  2003,  Oct  27.9972.  As  shown 
on  the  left  panel  of  Figure  9,  the  depth,  duration  and 
timing  of  this  event  were  consistent  with  an  egress  from 
transit  of  G1  876  c,  in  agreement  with  our  predictions. 
The  possibility  that  a  transit  had  been  observed  spurred 
us  to  orchestrate  a  distributed  observational  campaign 
during  the  2004  season. 

During  the  2004  June  24  opportunity,  observations 
were  attempted  by  ten  observers  in  Australia.  Poor 
weather  thwarted  all  but  one  observer,  and  no  transit  sig¬ 
nal  was  seen  in  the  data,  which  covered  about  half  of  the 
1  —  a  transit  window.  Eight  Australian  and  Japanese  ob¬ 
servers  participated  in  the  following  2004  July  24  planet 
c  campaign.  They  were  all  clouded  out.  Weather  im¬ 
proved  for  the  2004  August  23  planet  c  opportunity,  per¬ 
mitting  observations  by  five  of  the  nine  Australian,  South 
African,  and  German  observers.  No  transit  signal  was 
seen  with  photometric  depth  greater  than  1%. 

In  2004  October,  a  number  of  AAVSO  and  Transit- 
Search  observers  obtained  photometry  during  a  closely 
spaced  pair  of  transit  windows  for  planets  b  and  c.  Begin¬ 
ning  October  20  11:22  UT  and  ending  October  23  15:27 
UT,  a  total  of  2,795  CCD  observations  were  obtained.  In 
addition,  2,981  photometric  observations  were  obtained 
in  the  days  immediately  before  and  after  the  window  to 
set  baseline  activity  and  to  look  for  red  dwarf  flaring. 
Such  a  large,  distributed  network  of  observers  allowed 
coverage  of  the  window  with  very  few  observing  gaps  over 
the  3-day  period.  Gaps  in  the  coverage  were  smaller  than 
the  estimated  transit  period  or  forecast  window,  allowing 
us  to  conclude  that  transits  for  planets  b  and  c  did  not 
occur  during  the  October  2004  campaign.  Photometric 
data  from  the  campaign  are  plotted  in  Figures  9  and  10. 

Individual  observations  and  uncertainty  estimates  from 
the  October  2004  campaign  are  available  from  the 
AAVSO  web  site,  while  the  other  datasets  can  be  found 
at  the  TransitSearch  website14.  We  are  continuing  to 
obtain  photometry  of  G1  876,  with  the  goal  of  gaining 
a  better  understanding  of  the  photometric  variability  of 
mature  (age  exceeding  lGyr)  red  dwarfs. 

5.  DISCUSSION 

The  TransitSearch  photometric  network  has  been  in 
operation  for  three  years,  and  has  participating  observers 
capable  of  providing  fully  global  longitude  and  latitude 
coverage.  The  network  maintains  a  continuously  updated 
catalog15  of  the  known  census  of  extrasolar  planets.  The 
catalog  is  presently  the  only  available  source  for  (1)  a- 
priori  transit  probabilities  and  predicted  transit  depths, 
for  the  known  extrasolar  planets,  (2)  predicted  transit 
ephemeris  and  transit  window  uncertainties  based  on  or¬ 
bital  fits  to  published  radial  velocities,  and  (3)  the  results 
of  known  photometric  searches  for  transits  of  planet¬ 
bearing  stars.  To  date,  the  majority  of  negative  transit 
searches  have  not  been  reported  in  the  literature.  With 
this  paper,  we  hope  to  spur  a  reversal  of  that  trend. 

14  our  AAVSO  data  available  at 

http:/ /www. aavso.org/data/download/  while 

other  collected  data  can  be  found  at 
http:  / /www.  ucolick.org/  %7elaugh/G1876/ c.results.html 

10  available  at  http://www.ucolick.org/%7elaugh/ 


Due  to  significant  uncertainties  in  the  orbital  fits,  it  is 
generally  very  difficult  to  definitively  rule  out  transits  by 
planets  with  P  >  lOd.  The  G1  876  planets  form  a  no¬ 
table  exception  to  this  rule  because  of  their  high  signal- 
to-noise,  T  ~  50,  and  the  extensive  radial  velocity  data 
set  allows  for  a  very  accurate  orbital  fit.  Campaigns  by 
the  TransitSearch  collaboration  have  thus  far  led  to  tran¬ 
sits  of  a  number  of  other  shorter-period  planet-bearing 
stars  being  ruled  out  to  a  high  degree  of  confidence.  The 
planets  for  which  the  possibility  of  transits  have  been 
significantly  discounted  by  our  campaigns  include  HD 
217107  “b”  (P  =  7.127d,  Ptr  =  8.0%,  see  also  Vogt  et 
al.  2005),  HD  168746  “b”  (P  =  6.403d,  Ptr  =  8.1%,  see 
also  Pepe  et  al.  2002),  and  HD  68988  “b”  (P  =  6.276d, 
Ptr  =  9.0%,  see  also  Vogt  et  al.  2002). 

Furthermore,  planets  for  which  usable  photometry  of 
the  parent  stars  has  been  obtained  by  our  network  during 
the  3  —  (7  transit  windows,  but  for  which  transits  can 
not  yet  be  fully  ruled  out,  include  HD  188753Ab  (P  = 
3.348d,  PtI  =  11.8%),  HD  76700  “b”  (P  =  3.971d,  Ptr  = 
10.0%),  HD  13445  b  (P  =  15.76d,  Ptr  =  3.4%),  HD  74156 
“b”  (P  =  51.64d,  Ptr  =  3.8%),  HD  37605  “b”  (P  = 
54.23d,  Ptr  =  2.2%),  and  HD  80606  “b”  (P  =  111.4d, 
Ptr  =  1.7%).  Planets  for  which  campaigns  are  scheduled 
to  begin  soon  include  GL  581  “b”  (P  =  5.366d,  Ptr  = 
3.6%),  and  HD  99492  “b”  (P  =  17.04d,  Ptr  =  3.4%).  GL 
581  is  of  particular  interest  because  the  primary  star  is  an 
M  dwarf,  and  the  planet  (if  it  transits)  will  be  a  Neptune- 
mass  object  M sin(i)  =  0.052Mjup  (Bonfils  et  al.  2005). 
The  1.6%  transit  depth  would  allow  an  unprecedented 
physical  characterization  of  a  low-mass  extrasolar  planet. 
To  date,  our  G1  876  campaign  has  employed  the  most 
globally  far-reaching  distribution,  and  along  with  recent 
leaps  in  fidelity  in  similar  campaigns,  clearly  underscores 
the  utility  of  the  network. 

What  is  more,  G1  876  b  and  c  constitute  a  dramatic 
exception  to  the  emerging  aphorism  that  Jupiter-mass 
planets  are  rarely  associated  with  M  dwarfs.  Indeed, 
a  paucity  of  giant  planets  orbiting  red  dwarfs  seems  to 
be  a  natural  consequence  of  the  core  accretion  theory 
of  planet  formation  (Laughlin,  Bodenheimer  &  Adams 
2004).  The  core-accretion  theory  predicts,  however,  that 
Neptune-mass  and  smaller  mass  planets  will  be  very  com¬ 
mon  around  red  dwarfs.  This  stands  in  contrast  to  the 
predictions  of  Boss  (2006),  who  suggests  that  the  gravi¬ 
tational  instability  mechanism  is  the  dominant  mode  of 
giant  planet  formation.  If  this  is  true,  then  Jupiter-mass 
planets  should  be  just  as  common  around  red  dwarfs  as 
they  are  around  solar-type  stars.  Efforts  to  detect  tran¬ 
sits  of  small  planets  orbiting  small  stars  are  therefore 
likely  to  gain  increasing  importance  and  focus  over  the 
next  several  years. 

There  are,  in  fact,  several  observational  indications 
that  low-mass  planets  may  be  common  around  red 
dwarfs.  The  radial  velocity  surveys  have  recently  re¬ 
ported  the  detection  of  Neptune-mass  companions  in 
short-period  orbits  around  the  red  dwarfs  GL  581  and 
GL  436.  The  OGLE  team  of  observers  using  the  rni- 
crolensing  method  have  detected  the  signature  of  what 
appears  to  be  a  planet  with  5.5  Earth  Masses  orbiting 
a  distant  red  dwarf  (Bennet  et  al.  2006).  With  these 
points  in  mind,  we  can  put  the  results  of  our  G1  876 
campaign  in  a  broader  context.  A  ground-based  (and 
“fiscally  viable”)  photometric  network  such  as  Transit- 
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search,  especially  in  collaboration  with  a  dedicated  bank 
of  telescopes,  can  capably  monitor  individual  M  dwarfs 
to  search  for  transits  of  terrestrial-sized  planets.  Owing 
to  the  intrinsic  long-term  variability  of  M  dwarfs,  candi¬ 
date  planets  would  have  to  be  identified  on  the  basis  of 
a  full  transit  signature  in  a  single  time-series. 

A  0.1  Mq  M  dwarf  has  A*/i?0  ~  0.1,  Teff  ~  2750K, 
and  L*/Lq  ~  5  x  10~5;  a  habitable  planet  that  receives 
an  Earth-equivalent  flux  from  the  star  thus  needs  to  or¬ 
bit  at  a  distance  of  just  0.022  AU,  which  corresponds  to 
an  orbital  period  of  3.85  days.  An  Earth-like  planet  with 
such  a  period  would  be  rotationally  synchronized  to  the 
red  dwarf,  and  in  the  absence  of  any  significant  perturb¬ 
ing  bodies,  its  orbit  would  be  almost  perfectly  circular. 
Interestingly  however,  simulations  show  that  habitability 
is  unlikely  to  be  adversely  affected  by  a  spin-orbit  period 
synchronization.  Joshi,  Haberle,  &  Reynolds  (1997)  used 
a  global  climate  model  to  investigate  how  the  Earth’s  cli¬ 
mate  would  respond  if  the  Earth  were  tidally  locked  to 
the  Sun,  and  found  that  Earth  remains  habitable  in  this 
configuration  -  at  1  AU,  and  perhaps  throughout  the 
habitable  zone. 

The  M  =  7.5 M®,  P  =  1.9379  d  companion  to  G1 
876  demonstrates  that  it  is  not  unreasonable  to  expect 
terrestrial-mass  bodies  on  short-period  orbits  of  O.1M0 
stars.  Montgomery  &  Laughlin  (2006)  have  carried 
out  accretion  simulations  with  the  Wetherill-Chambers 
method  (Wetherill  1996,  Chambers  2001)  that  model  the 
accretion  of  terrestrial-mass  planets  in  short-period  or¬ 
bits  about  M  dwarfs.  These  simulations  provide  further 
evidence  that  terrestrial-mass  planets  will  commonly 
form  in  the  desired  short-period  orbits. 

If  habitable  planets  do  commonly  form  in  orbit  around 
low  mass  M  dwarfs,  the  chances  of  detecting  and  char¬ 
acterizing  them  are  surprisingly  good.  The  transit  of  an 
Earth-sized  planet  will  block  about  1%  of  the  stellar  flux 
of  a  0.1  M0  star.  For  a  planet  on  a  habitable  3.85  day 
orbit,  the  transit  will  be  relatively  brief,  ~  40  minutes. 
The  a-priori  geometric  probability  of  observing  such  a 
transit  is  2%. 

A  1%  photometric  dip  is  readily  detectable.  Amateur 
astronomers  who  participate  in  the  Transitsearch  collab¬ 
oration  routinely  achieve  detection  thresholds  of  consid¬ 
erably  better  than  1%,  as  evidenced  by  confirming  de¬ 
tections  of  HD  149026b  which  has  a  transit  depth  of  just 
0.3%.  Indeed,  capable  amateur  observers  have  demon¬ 
strated  the  photometric  capability  to  detect  the  passage 
of  a  Mars-sized  body  in  front  of  an  11th  magnitude  0.1 
solar  mass  red  dwarf. 

Inevitably  one  must  ask  how  many  suitable  red  dwarfs 
are  available  on  the  sky  as  a  whole.  Even  though  the 
lowest-mass  M  stars  are  the  most  common  type  of  star, 
they  are  also  exceedingly  dim.  A  dedicated  3m-class  tele¬ 
scope  would  be  required  to  properly  search  the  magni¬ 
tude  range  V=16-18,  where  observations  of  the  lowest- 
mass  red  dwarfs  begin  to  be  plentiful  (and  indeed,  many 
0.1  M0  stars  within  10  parsecs  still  likely  remain  to 
be  discovered).  In  the  meantime,  however,  we  recom¬ 
mend  that  observers  obtain  high-cadence  photometry  of 
nearby  single  stars  listed  in  the  RECONS  catalog16  of 
the  100  nearest  stellar  systems;  these  include  Proxima 

16  maintained  by  T.  Henry  and  collaborators  at  Georgia  State 
University. 


Centauri,  Barnard’s  Star,  Wolf  359,  Ross  154,  Ross  128, 
DX  Cancri,  GJ  1061,  GJ  54.1,  and  GJ  83.1.  We  esti¬ 
mate  that  each  one  has  a  ~  1%  chance  of  harboring  a 
detectible  transiting,  potentially  habitable  planet  -  read¬ 
ily  detectible  by  distributed  photometric  observation. 

6.  CONCLUSIONS 

Our  study  of  G1  876  has  produced  several  interesting 
results.  We  have  demonstrated  that  transits  of  planets 
b  and  c  are  not  currently  occurring,  and  we  have  thus 
modestly  constrained  the  parameter  space  available  to 
observationally  allowed  orbital  configurations  of  the  G1 
876  system.  The  absence  of  transits  is  consistent  with 
the  dynamical  study  of  Rivera  et  al.  (2005)  who  find 
that  a  co-planar  system  inclination  i  ~  50°  yields  the 
best  fit  to  the  radial  velocity  data.  Our  dynamical  anal¬ 
ysis  indicates  that  when  planetary  transits  in  a  strongly 
interacting  system  similar  to  G1  876  are  eventually  ob¬ 
served,  then  a  tremendous  amount  of  information  will  be 
obtained.  We  look  forward  to  the  discovery  of  the  first 
such  system. 

We  have  also  shown  that  a  globally  distributed  net¬ 
work  of  small-telescope  observers  can  effectively  provide 
definitive  photometric  follow-up  of  known  planet-bearing 
low-mass  stars.  Our  experience  with  G1  876  indicates 
that  intensive  photometric  monitoring  campaigns  of  in¬ 
dividual,  nearby  M  dwarfs  constitute  a  tractable  strategy 
for  planet  detection.  The  Sun’s  nearest  low-mass  stellar 
neighbors  should  therefore  be  scrutinized  very  carefully 
for  transits  by  short-period  terrestrial  planets.  A  low- 
budget  survey  of  this  type  provides  a  small,  but  never¬ 
theless  non-negligible  chance  to  discover  a  truly  habitable 
world  from  the  ground,  and  within  a  year  or  two.  Such  a 
detection  would  be  a  remarkable  and  exciting  discovery. 
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Fig.  1. —  Uncertainty  in  the  central  transit  time  for  planets  “b”  solid  line  and  “c”  dashed  line.  The  vertical  dashed  lines  delineate  the 
time  span  of  the  Keck  RV  observations  used  to  generate  the  co-planar  fits  that  form  the  basis  of  our  model  photometric  and  RV  curves. 
The  vertical  dotted  line  indicates  the  epoch  for  the  fits  (JD  2452490.0). 
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Fig.  2. —  Model  limb-darkened  light  curves  during  a  predicted  transit  by  planet  “b”  when  it  is  near  periastron  ( solid  curves)  and  near 
apastron  ( dashed  curves).  The  curves  correspond  to  different  assumed  inclinations  ( i  =  90°,  89.9°,  89.8°,  89.5°).  The  shallowest  occur 

for  i  =  89.5°. 
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Fig.  3. —  Model  transit  light  curves  for  planet  “c”,  when  it  is  near  periastron  ( solid  lines),  and  near  apastron  ( dashed  lines).  Inclinations, 
i  =  90°  through  i  =  89.0°  are  shown  in  decrements  of  0.1°.  Planet  “c’”s  eccentricity  e  ~  0.22  leads  to  a  greater  difference  in  transit  depth 
and  duration  between  periastron  and  apastron. 
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Fig.  4. —  Three-planet  RV  model  ( dashed  curve)  with  a  superimpo; 
planet  “c”  near  epoch  JD  2453301.79.  The  observed  Keck  RV  data  is  s' 
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Fig.  7. —  Model  i  =  89.3°  light  curves  for  G1  876  “b”  and  “c”  superimposed  on  Fairborn  Observatory  photometry.  Note  that  the 
intrinsic  long-term  variability  for  G1  876  over  a  ~  1000  day  period  is  of  order  0.05  magnitudes. 
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Fig.  8. —  Model  light  curves  for  G1  876  “c”  superimposed  on  Fairborn  Observatory  photometry  taken  near  three  individual  predicted 
transit  windows.  The  solid  lines  correspond  to  model  light  curves  arising  from  assumed  inclinations  running  from  i  =  90.0°  to  i  =  89.0°. 
The  light  lines  correspond  to  %  =  90°  model  light  curves  arising  from  bootstrap  fits  in  which  the  transit  is  3 a  early,  la  early,  and  lcr  and 
3a  late. 
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Fig.  9. —  Left:  Shankland  2003  data  for  “c”,  this  plate  shows  dynamical  model  light  curves  again  for  the  same  decrements  of  i ,  taken 
near  JD  2452940.4381;  relative  photometry,  which  began  in  astronomical  twilight  ( filled-in  circles)  Right:  Dynamical  model  light  curves 
for  G1  876  due  to  planet  “c”  for  i  =  90.0°  through  89.0°  in  0.1°  decrements,  taken  near  a  predicted  transit  centered  at  JD  2453301.7879. 
The  data  were  compiled,  baselined  and  shifted  to  employ  a  mean  magnitude  of  V=10.19,  based  on  relative  photometry  obtained  by  AAVSO 
and  TransitSearch  ( filled-in  circles).  The  light  lines  correspond  to  i  =  90°  model  light  curves  arising  from  bootstrap  fits  in  which  the 
transit  is  3cr  early,  lcr  early,  and  1  a  late  and  3 cr  late. 
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Fig.  10. —  Dynamical  model  light  curves  for  G1  876  due  to  planet  “b”  for  i  =  90.0°  through  89.5°  in  0.1°  decrements,  near  a  predicted 
transit  centered  at  JD  2453300.3806.  As  with  the  previous  figure,  the  data  were  baselined  and  shifted  to  employ  a  mean  magnitude  of 
V=10.19,  based  on  relative  photometry  obtained  by  AAVSO  and  Transits  earch(filled-in  circles).  The  light  lines  correspond  to  i  =  90° 
model  light  curves  arising  from  bootstrap  fits  in  which  the  transit  is  3cr  early,  la  early,  and  la  late  and  3 a  late. 
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TABLE  1 

Co-Planar  Three-Planet  i  =  90°  Jacobi  Parameters  of  the  Gl  876  System 


